BELLE-CONF-1271 



Evidence for a Z^(10610) in Dalitz analysis of 
Y(5S') r{nS)7TV 

I. Adachi/^ K. Adamczyk,'^^ H. Aihara,^'^ K. Arinstein,^ Y. Arita,^^ D. M. Asner,^^ 
T. Aso/^ V. Aulchenko,^ T. Aushev,^^ T. Aziz,^^ A. M. Bakich,^^ y. Ban,^^ Barberio,^^^ 
M. Barrett/2 A. Bay,^^ 1. Bedny,^ M. Belhorn,^ K. Belous,^^ V. Bhardwaj,^^ B. Bhuyan/^ 
M. Bischofberger,^^ S. Blyth,^^ A. Bondar,^ G. Bonvicini/^ A. Bozek,^^ M. Bracko,^^'^^ 
J. Brodzicka,^^ O. Brovchenko,^^ T. E. Browder/^ M.-C. Chang,'' P. Chang,^^ Y. Chao,^^ 

V. Chekelian,36 A. Chen/2 K.-F. Chen,^^ P. Chen,^^ B. G. Cheon/^ K. Chilikm,^^ 
R. Chistov,24 I.-S. Cho/^ K. Chc^^ K.-S. Choi/^ S.-K. Choi/° Y. Choi,*'^ j Crnkovic/^ 

J. Dalseno,36'65 m. Danilov,^^ J. Dingfelder/ Z. Dolezal,^ Z. Drasal,^ A. Drutskoy,^^ 
W. Dungel,20 D. Dutta/^ S. Eidelman,^ D. Epifanov,^ S. Esen,^ J. E. Fast,^! M. Feindt,^^ 
M. Fujikawa,^^ V. Gaur,^^ N. Gabyshev,^ A. Garmash,^ Y. M. Goh," B. Golob,33'25 
M. Grosse Perdekamp/^'^^ H. Guo,^^ J. Haba/^ y. L. Han/^ K. Hara/^ T. Hara/^ 
Y. Hasegawa,^^ K. Hayasaka,^^ H. Hayashii,^^ D. Heffernan,^° T. Higuchi/^ Y. Horii,'^^ 
Y. Hoshi,^^ K. Hoshina/^ W.-S. Hou,^^ Y. B. Hsiung,^^ H. J. Hyun,^! Y. Igarashi/^ 
T. lijima,^^'^^ M. Imamura,^* K. Inami,^^ A. Ishikawa,^^ R. Itoh/^ M. Iwabuchi/'^ 
M. Iwasaki,^^ Y. Iwasaki/^ T. Iwashita,''^ S. IwataJ^ I. Jaegle/^ M. Jones/^ T. Julius,^^ 
D. H. Kah,3i H. KakunoJ^ J. H. Kang/^ P. Kapusta,^^ S. U. Kataoka,^° N. Katayama/^ 
H. Kawai,^ T. Kawasaki,^^ H. Kichimi/^ C. Kieslmg,^^ B. H. Kim,^^ H. J. Kim,^^ 
H. O. Kim,3i J. B. Kim,^^ J. H. Kim,^^ K. T. Kim,^^ M. J. Kim,^^^ S. K. Kim,^^ 
Y. J. Kim,28 K. Kinoshita,^ J. Klucar,^^ B. R. Ko,^^ N. Kobayashi/o S. Koblitz,^^ 
P. Kodys,3 Y. Koga,38 S. Korpar,^^^' R. T. Kouzes,^^ M. Kreps,^^ P. Krizan,33.25 
P. Krokovny,^ B. Kronenbitter,^'' T. Kuhr,^'' R. Kumar, T. Kumita,^^ E. Kurihara,^ 
Y. Kuroki,^° A. Kuzmin,^ P. Kvasnicka,^ Y.-J. Kwon,^'' S.-H. Kyeong,''^ J. S. Lange,^ 
M. J. Lee,6° S.-H. Lee,^^ M. Leitgab,^^'^^ R .Leitner,^ J. Li,^^ X. Li,^^ Y. LiJ^ J. Libby,^^ 
C.-L. Lim,^^ A. Limosani,^^ C. Liu,^^ Y. Liu,^ Z. Q. Liu,^^ D. Liventsev,^^ R. Louvot,^^^ 
J. MacNaughton,i3 j3_ Marlow,^"^ D. Matvienko,^ A. Matyja,"^^ S. McOnie,^^ Y. Mikami,^^ 
K. Miyabayashi,^! Y. Miyachi,™ H. Miyata,^^ Y. Miyazaki,^^ R. Mizuk,^^ G. B. Mohanty,''^ 

D. Mohapatra,^^ A. Moll,^^'^^ T. Mori,^^ T. Miiller,^^ N. Muramatsu,^^ R. Mussa,^^ 
T. Nagamine,^^ Y. Nagasaka,^'^ Y. Nakahama,^^ I. Nakamura,^^ E. Nakano,'*^ T. Nakano,^^ 
M. Nakao,^^ H. Nakayama,^^ H. Nakazawa,''^ Z. Natkaniec,^^ M. Nayak,^^ E. Nedelkovska,^^ 
K. Negishi,68 K. Neichi,^^ S. Neubauer,^^ C. Ng,^^ Niiyama,=^° S. Nishida,^^ 
K. Nishimura,i2 q. Nitoh,^^ Nozaki,^^ Qgawa,^^ S. Ogawa,^^ T. Ohshima,^^ 
S. Okuno,26 S. L. Olsen,^^'!^ Y. Onuki,^^ W. Ostrowicz,^^ H. Ozaki,!^ R Pakhlov,^^ 
G. Pakhlova,24 H. Palka,^^ C. W. Park,^^ Park,^! H. K. Park,^! K. S. Park,^^ 
L. S. Peak,63 K. Pedlar,^^ T. Peng,^^ R. Pestotnik,^^ M. Peters,^^ m. Petric,^^ 
L. E. Piilonen,'^^ A. Poluektov,^ M. Prim,^^ K. Prothmami,^^'^^ B. Reisert,^*' M. Ritter,^^ 
M. R5hrken,27 J. Rorie,^^ m. Rozanska,^^ S. Ryu,^° H. Sahoo,^^ k. Sakai,!^ Y. Sakai,^^ 
S. Sandilya,^^ D. Santel,^ L. Santelj,^^ T. Sanuki,^^ N. Sasao,30 Y. Sato,''^ O. Schneider,^^ 
P. Schonmeier,^^ C. Schwanda,^° A. J. Schwartz,^ R. Seidl,^^ A. Sekiya,^^ K. Senyo,''^ 
O. Seon,38 M. E. Sevior,^^ L. Shang,^^ M. Shapkm,^! V. Shebalin,^ C. P. Shen,^^ 
T.-A. Shibata,^° H. Shibuya,^^ S. Shinomiya,^^ J.-G. Shiu,^^ B. Shwartz,^ A. Sibidanov,^^ 
F. Simon,36.65 j. b. Singh,^^ R. Smha,^^ P. Smerkol,^^ Y.-S. Sohn,^^ A. Sokolov,^! 



E. Solovieva,^^ S. Stanic,^^ M. Staric,^''' J. Stypula,^''^ S. Sugihara,^^ A. Sugiyama,"''^ 
M. Sumihama,^ K. Sumisawa/'^ T. Sumiyoshi/-'^ K. Suzuki, '^^ S. Suzuki, S. Y. Suzuki, 
H. Takeichi,^^ U. Tamponi,^'^ M. Tanaka,^"^ S. Tanaka,^^ K. Tanida,^° N. Taniguchi,^^ 

G. Tatishvili,^^ G. N. Taylor, Y. Teramoto,'*^ F. Thome, I. Tikhomirov,^^ K. Trabelsi,^^ 

Y. F. Tse,^^ T. Tsuboyama,^^^ M. Uchida,^° T. Uchida,^^ Y. Uchida,'^ S. Uehara,^^ 
K. Ueno,^^ T. Uglov,^^ Y. Unno,!^ S. Uno,^^^ P. Urquijo,^ Y. Ushiroda,^^ y. Usov,^ 
S. E. Vahsen,^^ P. Vanhoefer,^^ G. Varner,^^ K. E. Varvell,^^ K. Vervink,^^ A. Vinokurova,^ 

V. Vorobyev,2 A. Vosscn,i« C. H. Wang,^''^ J. Wang,^^ ^ .2. Wang,^^ P. Wang,!^ 
X. L. Wang,i9 M. Watanabc,^^ Y. Watanabe,^^ R. Wedd,^^ E. White,^ J. Wicht,^^ 

L. Widhalm,20 J. Wiechczynski,^^ K. M. Williams,^^ E. Won,^^ B. D. Yabsley,^^ 

H. Yamamoto,^^ J. Yamaoka,^^ Y. Yamashita,^^ M. Yamauchi,^^ C. Z. Yuan,^'^ Y. Yusa,^'^ 
D. Zander,27 C. C. Zhang,^^ L. M. Zhang,^^ Z. P. Zhang,^^ L. Zhao,^'^ V. Zhilich,^ 

P. Zhou,^^ V. Zhulanov,^ T. Zivko,^^ A. Zupanc,^^ N. Zwahlen,^^ q Zyukova^ 

(The Belle Collaboration) 

^ University of Bonn, Bonn 
^Budker Institute of Nuclear Physics SB RAS and 
Novosibirsk State University, Novosibirsk 630090 
^Faculty of Mathematics and Physics, Charles University, Prague 
^Chiba University, Chiba 
^University of Cincinnati, Cincinnati, Ohio 45221 
^Department of Physics, Fu Jen Catholic University, Taipei 
Justus-Liebig-Universitdt Ciefien, Giefien 
^Gifu University, Gifu 
^The Graduate University for Advanced Studies, Hayama 
^^Gyeongsang National University, Chinju 
^^Hanyang University, Seoul 
University of Hawaii, Honolulu, Hawaii 96822 
^^High Energy Accelerator Research Organization (KEK), Tsukuba 

^■^ Hiroshima Institute of Technology, Hiroshima 
University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 
^^Indian Institute of Technology Cuwahati, Cuwahati 
^''Indian Institute of Technology Madras, Madras 
Indiana University, Bloomington, Indiana 4.74.08 
^^Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
^^Institute of High Energy Physics, Vienna 
Institute of High Energy Physics, Protvino 
Institute of Mathematical Sciences, Chennai 
^^INFN - Sezione di Torino, Torino 
^'^Institute for Theoretical and Experimental Physics, Moscow 
^^J. Stefan Institute, Ljubljana 
Kanagawa University, Yokohama 
^'^Institut fiir Experimentelle Kernphysik, 
Karlsruher Institut fiir Technologic, Karlsruhe 
Korea Institute of Science and Technology Information, Daejeon 



Korea University, Seoul 
^'^Kyoto University, Kyoto 
Kyungpook National University, Taegu 
^^Ecole Polytechnique Federals de Lausanne (EPFL), Lausanne 
Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 
^'^ Luther College, Decorah, Iowa 52101 
University of Maribor, Maribor 
^'^ Max- Planck- Institut fiir Physik, Miinchen 
^''^ University of Melbourne, School of Physics, Victoria 3010 
Graduate School of Science, Nagoya University, Nagoya 
Kobayashi-Maskawa Institute, Nagoya University, Nagoya 
^^Nara University of Education, Nara 

^^Nara Women's University, Nara 
'''^National Central University, Chung-li 
^'^ National United University, Miao Li 
"^"^ Department of Physics, National Taiwan University, Taipei 
^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 
^^Nippon Dental University, Niigata 
^''Niigata University, Niigata 
'^^ University of Nova Gorica, Nova Gorica 
Osaka City University, Osaka 
Osaka University, Osaka 
^^Pacific Northwest National Laboratory, Richland, Washington 99352 
^^Panjab University, Chandigarh 
Peking University, Beijing 
Princeton University, Princeton, New Jersey 08544 
Research Center for Electron Photon Science, Tohoku University, Sendai 
Research Center for Nuclear Physics, Osaka University, Osaka 
^^RIKEN BNL Research Center, Upton, New York 11973 
^^Saga University, Saga 
University of Science and Technology of China, Hefei 
^'^Seoul National University, Seoul 
^^Shinshu University, Nagano 
Sungkyunkwan University, Suwon 
^^School of Physics, University of Sydney, NSW 2006 
Tata Institute of Fundamental Research, Mumbai 
^^Excellence Cluster Universe, Technische Universitdt Miinchen, Garching 

^^Toho University, Funabashi 
^'^Tohoku Gakuin University, Tagajo 
Tohoku University, Sendai 
^^Department of Physics, University of Tokyo, Tokyo 
™ Tokyo Institute of Technology, Tokyo 
'^^ Tokyo Metropolitan University, Tokyo 
'''^ Tokyo University of Agriculture and Technology, Tokyo 
''^Toyama National College of Maritime Technology, Toyama 
'^^CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
'^^ Wayne State University, Detroit, Michigan 48202 



Yamagata University, Yamagata 
^^Yonsei University, Seoul 

Abstract 

We report the first observation of T(5S') — t- T(1, 2S')7r''7r*' decays. Evidence for the Z^(10610) 
with 4.9(7 significance is found in a Dalitz plot analysis of T(5S') — )• T(25)7r'^7r'' decays. The results 
are obtained with a 121.4fb~^ data sample collected with the Belle detector at the T(5S') resonance 
at the KEKB asymmetric-energy e"^e~ collider. 



PACS numbers: 14.40.Pq, 13.25. Gv, 12.39.Pn 



INTRODUCTION 



Recently the Belle Collaboration reported the observation of two narrow structures in 
7r^T(?T,S') invariant mass in the T{5S) — )■ T(nS')7r+7r~ decays (n = 1, 2, 3) The measured 
masses and widths of the two structures are Mi = 10607.2±2.0 MeV/c^, Ti = 18.4±2.4 MeV 
and Ma = 10652.2 ± 1.5 MeV/c^ = 11.5 ± 2.2 MeV, respectivelv. Angular analysis sug- 
gests that these states have I^{J^) = 1^(1^) quantum numbers j2|. The measured masses 
are a few MeV/c^ above the thresholds for the open beauty channels B*B (10604.6 MeV/c^) 
and B*B (10604.6 MeV/c^) suggesting a "molecular" nature for these states, which is con- 
sistent with many of their observed properties 0. This observation motivates us to search 
for a neutral partner of these states in the resonant substructure of T(5S') — ?■ T(r?,5')7r°7r° 
decays. 

SELECTION CRITERIA 

We use a 121.4 ± 1.9 fb^^ data sample collected on the peak of the T(5S') resonance with 
the Belle detector j3] at the KEKB asymmetric energy e~^e~ collider (Ef. The Belle detector is 
a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector, a central 
drift chamber, an array of aerogel threshold Cherenkov counters, a barrel-like arrangement of 
time-of-flight scintillation counters, and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented to detect mesons and to identify 
muons. The detector is described in detail elsewhere 

T(5S') candidates are formed from T(nS')7r°7r'' (n = 1,2) combination. We recon- 
struct T{nS) candidates from pairs of leptons (e^e~ and fJ'~^fi~) with invariant mass in 
the range from 8 to 11 GeV/c^. An additional decay channel is used for the T(2S'): 
T(2S') — )■ T(lS')[/"''Z~]7r"*'7r~. Charged tracks are required to have transverse momentum, 
Pt >, greater than 50MeV/c. We also impose a requirement on the impact parameters: 
dr < 0.3 cm and \dz\ < 2.0 cm, where dr and dz are the impact parameters in the r- 
(j) and longitudinal directions, respectively. Muon and electron candidates are required to 
be positively identified. No requirement on the particle identification is used for the pi- 
ous. Candidate vr" mesons are selected from pairs of photons with an invariant mass within 
15 MeV/c^ of the nominal 7r° mass. Energy greater than 50 (75) MeV is required for each 
photon in the barrel (endcap). We use the quality of the vr" mass-constrained fits to sup- 
press the background; the sum of x^(^i) + X^(^2) is required to be less than 20 (10) for the 
T{nS) T(l^)7r+7r- {T{nS) -> e+e-). 

We use the energy difference, AE = -Ecand — Ecu, and momentum P to suppress back- 
ground, where -Ecand and P are the energy and momentum of the reconstructed T(5S') 
candidate in the center-of-mass (cm.) frame, and -Ecm is the cm. energy of the two 
beams. T(5S') candidates are required to satisfy the requirements —0.2 < AE < 0.14 GeV 
and P < 0.2 GeV/c. The large potential background from QED processes such as 
e~^e^ —7- l^l^{n)'y is suppressed using the missing mass associated with the system, 
Mmiss(^^^"), calculated as Mmiss (/+/") = {Ecm - Ei+i-Y - P,^+,_, where Ei+i- and Pi+i- 
are the energy and momentum of the /"'"Z" system measured in the cm. frame. We require 
Mmiss(/+/-) > 0.15(0.30) GeV/c2 for the T{nS) fi+fi' (e+e"). We select the candidate 
with the smallest x^(^i) + X^(^2) i^ ^^e rare cases (1-2%) in which there is more than one 
candidate in the event. 
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FIG. 1: The tt^tt^ missing mass distribution for T(n5)7r'^7r'', (a) T(n5) — )• fJ^'^fJ.^ and (b) T{nS) — >■ 
e^e~ candidates. The M(T(lS')7r"'"7r~) distribution for T{2S) — )• T(15)7r+7r~ candidates is shown 
in (c). Histograms represent the data. The sohd curves show the fit result while the dashed curves 
correspond to the background contributions. 

T(55') — 7- T{nS)[l^l^]'n-^7T^ candidates are identified via the missing mass recoiling against 
the 7r°7r° system, Mmiss(7r°7r°). Figures[T](a) and[T](b) show the Mmiss(vr°7r°) distributions for 
T{5S) — i- T {nS)[l^l^]TT^iT^ candidates. We fit these distributions to extract the T{nS) signal 
yield. The signal probability density function (PDF) is described by a sum of two Gaussians 
for each T{nS) resonance with parameters fixed from the signal Monte Carlo (MC) sample. 
The background PDF is parameterized by the sum of constant and exponential functions. 

For T{5S) — i- T(2S')[T(lS')7r+7r^]7r°7r° decays, T(IS') candidates are selected from 
pairs with invariant mass within 150MeV/c^ of the nominal T(IS') mass. A mass- 
constrained fit is used for T(IS') candidates to improve the momentum resolution. We apply 
the same requirements on AE and P described above for reconstructed T(5S') candidate. 
We use the invariant mass of T{lS)7i^n^ to select the signal candidates. Figure [T](c) shows 
the M(T(15')7r+7r^) distribution for the [T(lS')7r"'"7r^]7r'^7r° events. We fit this distribution 
to extract the T(25')[T(lS')7r"'"7r~] signal yield. The signal PDF is described by a Gaussian 
function with parameters fixed from signal MC. The background PDF is described by a 
constant. The cross-feed from the decay T(55') — )■ T(2S')[T(15')7r°7r°]7r+7r^ contributes as a 
broad peak around 10.3 GeV/c^. Its contribution is parameterized by a Gaussian function. 

Table |T] summarizes the definition of the signal region, signal yield, MC efficiency, mea- 
sured branching fraction (only the statistical uncertainty is shown), number of selected events 
and purity. The reconstruction efficiency is obtained using MC with with the T(nS')7r°7r° sys- 
tem distributed uniformly over three-body phase space. The branching fraction is calculated 
by B = g£j(e+e-^T(55)) ' '^here iVgig is number of signal events, e is reconstruction efficiency, £ 
is integrated luminosity. We use the value of a{e~^e~ — )■ T{5S)) = 0.340 ±0.016 nb obtained 
with 121.4 fb"^ data. 

The main sources of systematic uncertainty in the branching fraction measurement are: 
uncertainties in the signal and background PDFs used in the fit mainly due to data/MC 
width differences: 5%; uncertainty in the e~^e~ — )■ T{5S) cross section: 5%; B(T{nS) — )■ 
/+/-): 2% and 9% for the T(l^) and T{2S) [7]; luminosity: 1.5%; tt" reconstruction: 5%; 
muon identification: 1%; electron identification: 3%; tracking: 0.7%. The total systematic 
errors are 9% for T(lS')7r*^7r° and 13% for T(2S')7r°7r°. We calculate the weighted average 
of B(T{5S) — i- T(nS')7r*^7r°) in the various T{nS) decay channels and obtain B{T{5S) 



TABLE I: Definition of the signal region, signal yield, MC efficiency, measured branching fraction, 
number of selected events and purity. 



Final state 


Signal region, GeV/c^ 


Signal yield 


e, % 


B, 10-3 


Events Purity 


T(15) - 




9.41 < MmissCvrOvrO) < 9.53 


261 ± 15 


11.2 


2.28 ±0.13 


247 


0.95 


T(15) - 




9.41 < Mmiss(^°vrO) < 9.53 


123 ± 13 


5.61 


2.15 ±0.23 


140 


0.78 


T{2S) - 


fi+fi- 


9.99 < Mmiss(7r07r°) < 10.07 


241 ± 18 


8.04 


3.77 ±0.28 


253 


0.87 


T{2S) - 


-> e+e~ 


9.99 < Mmiss(vrO^O) < 10.07 


108 ± 13 


3.58 


3.84 ± 0.46 


151 


0.66 


T{2S) - 


T(15)7r+^- 


10.00 < M(T7r+7r-) < 10.05 


24 ±5 


2.27 


2.85 ± 0.60 


28 


0.86 
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FIG. 2: Dalitz plots for selected (a) T(lS)7r'^7r'', (b) T(25)7r''7r'' candidates. Dalitz plots for events 
in the (c) T{lS)TT^7r^, (d) T{2S)tt^7:^ sidebands. 



T(15)7rV) = (2.25 ± 0.11 ± 0.20) x 10-^ and i3(T(55) ^ T(25)7r07r0) = (3.66 ± 0.22 ± 
0.48) X lO"'^. These are approximately one half of the corresponding values of B{T{5S) — )■ 
T{nS)7i^7i~) [9'], consistent with expectations from isospin. 



DALITZ ANALYSIS 

We define the following sideband regions for the study of background: for the T(15')7r'^7r° 
final state : 9.20GeV/c2 < M^i^^in'^'K^) < 9.35GeV/c2 and 9.60GeV/c2 < M^i^J^Tc^n^) < 
9.75GeV/c2; and for the T(25)7rV final state: 9.80GeV/c2 < M^i,,{7i^n^) < 9.95GeV/c2 
and 10.15 GeV/c^ < M„,iss(7r°7r°) < 10.30 GeV/c^. Figure |2] shows the Dalitz plot distribu- 
tions for the selected T(55') — T{nS)7r^7r^ candidates from the signal region and sidebands. 
Before analyzing Dalitz distributions for events in the signal region, we determine the PDF 
for background. Samples of background events are selected in T{nS) mass sidebands and 
then refitted to the nominal mass of the corresponding T{nS) state to match the phase 
space boundaries. We parameterize the background PDF by the following function: 

1 ±Piexp(-giS3 ±p2exp(-g2(smm - ^2)), (1) 

where pi, p2, qi and q2 are fit parameters. Here S3 = M^(7r°7r°) and Smin = Min(si,S2), 
= M2(T(n5)7r°2). The kinematical limit as is 92 and 103 GeVVc^ for the T(15) 
and T{2S), respectively. Variation of the reconstruction efficiency over the Dalitz plot is 
determined using a large sample of MC with a uniform phase space distribution. We use 
the following function to parameterize efficiency variations: 

e = 1 ± C{1 - e-("3-ao)/feo||l _ g-(ai-Sn,ax)/fel|^ (2) 



where c, bo and bi are fit parameters. Here Smax = Max(si,S2). The parameters Qq and ai 
are defined as Oq = 4m^05 c^i = {^^{53) ~ mj,oY. 

The amphtude analysis of three-body T(5S') — )■ T{nS)7i^n^ decays uses an unbinned 
maximum likehhood fit. We describe the three-body signal amplitude as a sum of quasi- 
two-body amphtudes: 

M(si, S2) = Azi + Az2 + Af^ + Af^ + A^, , (3) 

where Azi and Az2 are amplitudes for contributions from the Z^(10610) and Z^(10650), 
respectively. The amplitudes Af^, Af^ and A^,. are the contributions from the tc^tt^ system 
in an /o(980, /2(1275) and a non-resonant state, respectively. Here we assume that the 
dominant contributions are from amplitudes that preserve the orientation of the spin of 
the heavy quarkonium state and thus, both pions in the cascade decay T(5S') — )■ Z^n^ — )■ 
T{nS)n^n^ are emitted in an S'-wave with respect to the heavy quarkonium system. As 
demonstrated in Ref. [2|, angular analysis supports this assumption. Consequently, we 
parameterize both amplitudes with an iS-wave Breit-Wigner function 



BW(s,M,T) = , (4) 

^ ' ' ^ M2 - s - iMT ' ^ ^ 

where we neglect the possible s dependence of the resonance width. Both amplitudes are 
symmetrized with respect to 7r° interchange. The masses and widths are fixed to the values 
obtained in the T(n^)7r+7r- analysis: M{Zi) = 10607.2 MeV/c^, T{Zi) = 18.4MeV/c, 
M{Z2) = 10652.2 MeV/c^, 1(^2) = 11.5MeV/c [1]. Contributions from the /o(980) and 
/2(1275) are also included in the fit. We use a Flatte function for the /o(980) and a Breit- 
Wigner function for the /2(1275). Coupling constants of the /o(980) were fixed at values 
from the B+ K+n+n- analysis: M = 950MeV/c2, g^^ = 0.23, qkk = 0-73 i6]. The mass 
and width of the /2(1275) resonance are fixed to the world average values [7|. Following 
suggestions in Ref. j8|, the non- resonant amplitude A^r is parameterized as 

A^, = A^^e^^- + Ale'^^'^ss , (5) 

where Aj^., A^^., and (p^^ are free parameters in the fit. As there is only sensitivity to the 
relative amplitudes and phases between decay modes, we fix A^^^. = 10.0 and 0^^. = 0.0. 
The logarithmic likelihood function is defined as 

C = -2J2 log{e(si, S2)(/sig^(si, S2) + (1 - /sig)5(si, S2))} , (6) 

where S{si,S2) is |M(si,S2)P convoluted with the detector resolution (6.0MeV/c2 for 
T{nS)7r^ combinations), e(si,S2) describes the variation of the reconstruction efficiency 
over the Dalitz plot and fsig is the fraction of signal events in the data sample. The frac- 
tion /sig is determined separately for each T{nS) decay mode (see Table [T]). The function 
B{si,S2) describes the distribution of background events over the phase space. Both prod- 
ucts e(si, S2) ■ S{si, S2) and e(si, S2) ■ B{si, S2) are normalized to unity. 

Results from one-dimensional projections of the fits are shown in Figs. [3] and HI These 
projections are very similar to the corresponding distributions in T(nS')7r+7r^ [1]. A Z^ 
signal is most clearly seen in M(T7r'^)jnax- Table HTl shows the values and errors of amplitudes 
and phases obtained from the fit to the T(nS')7r°7r° Dalitz plot. The statistical significance 
of the Z°(10610) signal in the T(2S')7r°7r° sample is 5.3cr. This value is obtained from the 




FIG. 3: Comparison of tlie fit results (open histograms) witli experimental data (points with error 
bars) for T(lS')7r°7r° events in the signal region. Red and blue open histograms show the fit with 
and without Z^'s, respectively. Hatched histograms show the background components. 




FIG. 4: Comparison of the fit results (open histograms) with experimental data (points with error 
bars) for T(2S')7r°7r° events in the signal region. Red and blue open histograms show the fit with 
and without Z^'s, respectively. Hatched histograms show the background components. 



TABLE II: Results of the Dalitz plot fit of T(n5)7r°7r° events. 



Model 


T(15)7r07r0 
with Z°'s 


T(15)7r07r0 
w/o Z°'s 


with Z°'s 


T(25)7r07r0 
with Zy only 


T(25)7r07r0 
w/o Z°'s 


A(Z?) 
A{Zl) 


-36 ± 50 
0.60l°i^ 
-59 ± 60 


0.0 (fixed) 
0.0 (fixed) 


0.58l°:?i 

-113 ± 14 
n 07+0.20 

'-'••J' -0.16 

-125 ±27 


47+0-15 
-0.11 

-117 ± 17 
0.0 (fixed) 


O.O(fixed) 
0.0 (fixed) 


A{f2) 


15.7 ± 2.0 


14.6 ± 1.6 


18.2 ± 7.3 


23.9 ±7.3 


28.2 ± 7.0 


0(/2) 


60± 11 


51 ±9 


36 ±21 


28± 13 


28± 10 


Mh) 


1.07 ±0.15 


0.97 ±0.12 


11.5 ± 1.9 


10.5 ± 1.9 


8.2 ±2.1 


m) 


168 ± 11 


163 ± 10 


211 ±6 


213 ±7 


210 ±8 


^nr 


15.2 ± 1.2 


13.9 ±0.7 


34.7 ±4.9 


31.8 ±4.3 


24.6 ±4.2 




162 ±4 


161 ±4 


80± 12 


85± 13 


93± 15 


-21og£ 


-316.7 


-312.4 


-193.1 


-186.6 


-154.5 



TABLE III: Two solutions found in the Dalitz plot fit of T{2S)tt^tt^ events. 



Solutions 


w/o 
A 


w/o Z° 
B 


with Zl 
A 


with Zf 
B 


with Z°'s 
A 


with Z°'s 
B 


0(^1°) 


0.0 (fixed) 

— 

0.0 (fixed) 


0.0 (fixed) 

— 

0.0 (fixed) 


46+°-^^ 
-117 ± 14 
0.0 (fixed) 


1 qc+0.64 
"3-0.33 

88± 18 
0.0 (fixed) 


0.58l°:?l 
-113 ± 14 

n Q7+0.20 
^•"5' -0.16 

-125 ±27 


1.42 ±0.48 
91 ±21 

0.66 ± 0.40 
124 ± 37 


0(/2) 

Aifo) 

42 

nr 


OQ 1 ^ n 
ZO.Z ± /.U 

28± 10 
8.2 ±2.1 

210 ±8 
24.6 ± 4.2 

93± 15 


41. ± y.u 

-1 ± 14 

13.3 ±3.6 
131 ± 11 
44.2 ± 10.1 

-70 ± 16 


23. y ± 1.6 

18± 13 
10.5 ± 1.9 

213 ± 7 
31.8 ±4.3 

85± 13 


10± 16 
13.4 ±4.2 
134 ± 15 
50.4 ± 12.2 
-69 ± 22 


io.z ± /.3 
36± 21 

11.5 ± 1.9 
211 ±6 

34.7 ±4.9 
80± 12 


43.0 ± 15. D 

132 ± 19 
12.6 ±4.9 
132 ± 19 
50.8 ± 13.7 

-72 ± 25 


-21og£ 


-154.5 


-155.4 


-186.6 


-186.3 


-193.1 


-191.2 


p- value, A( 


-21og£), 


with two degrees 


. of freedom, 


i.e. A{-2\og 


;£) = 21og 


;£(w/oZO)- 



2 log £(10610). We also perform a fit with the Z°(10610) mass as a free parameter. The fit 
with both Z°'s gives M(Z0(10610)) = 10609^^ MeV/c^. A similar value, 10603 ± 6 MeV/c^, 
is obtained in a fit in which only the Z°(10610) is included. The signal for the Z^(10610) 
is not significant in the fit to the T(lS')7r°7r° events due to the smaller relative branching 
fraction. The signal for the Z^(10650) is not significant in either T{1, 23)71^7:^ dataset. Our 
data do not contradict the existence of Z°(10650), but the available statistics is not enough 
for the observation of this state. 

We search for multiple solutions by doing one thousand fits with randomly assigned 
amplitudes and phases taken from a model without contributions. We find an additional 
solution in the T{2S)7i^tt^ final state and no other solutions in T{lS)n^7r^. Table UTTl shows 
the values and errors of amplitudes and phases obtained for both solutions. The second 
solution is referred to "Solution B". The Z°(10610) significance in Solution B is 5.3cr, almost 
the same as in Solution A (the baseline fit). 

STUDY OF SYSTEMATIC UNCERTAINTIES IN DALITZ ANALYSIS 

Experimental errors may arise from the uncertainty in parameterization of the background 
PDF. We determine this uncertainty by varying parameters of the background PDF. We 
use different sideband sub-samples to determine PDF parameters: the low-mass sideband 
only, or the high-mass sideband, or T{nS) — )■ e~^e~ events only, or T{nS) — fi~^fj^~ events 
only. The statistical significance of the Z°(10610) in all fits is greater than 4.9o". Another 
source of systematic uncertainty is the determination of signal efficiency. To estimate this 
effect we perform two additional fits with a modified efficiency function: -\/e(s^^^^I7s3) and 
e^''^(smax, S3). The result, A(— 21og£) for models with Z°(10610) and without Z°'s changes 
from 32.1 to 32.9 and 31.2, respectively, the difference in Z^(10610) significance is less than 
0.1. We also perform a fit with a modified detector resolution function: the resolutions are 
varied from 4 to 8MeV/c^ instead of the nominal 6MeV/c^ to take into account the effect 



of different momentum resolutions in MC and data. The resulting changes in A(— 21og£) 
are less than 0.5 (difference in Z^(10610) significance is less than 0.05). 

The model uncertainty originates mainly due to the parameterization of the non-resonant 
amplitude. To estimate it we vary the model used to fit the data. Three additional models, 
based on solution A, are used: with an additional (t(600) resonance, parameterized by a 
Breit-Wigner function with M = GOOMeV/c^ and T = 400 MeV/c, a model with A^r = 
ae*"^" + be''l"'^/s(^l^7l^, and a model without any /o(980) contribution (to fit T(25')7r°7r° 
events only). The smallest Z^(10610) significance is obtained in the last model: 4.9(j. 
We use this value as the final Z°(10610) significance. Fits with the Z°(10610) mass as a 
free parameters give values from 10603 to 10615 MeV/c^. We use ±6MeV/c^ as a model 
uncertainty for the Z°(10610) mass. 

CONCLUSION 

We report the observation of T(5S') — t- T{1,2S)7i^7t^ decays. The measured branching 
fractions, B{T{5S) T(15)7rV) = (2.25 ± 0.11 ± 0.20) x 10"^, i3(T(5S) ^ T(2^)7rV) = 
(3.66 ±0.22 ±0.48) x 10"'^, are found to be consistent with the the expectation from isospin, 
scaling from i3(T(5^) ^ T(n^)7r+7r-) [91]. 

Evidence of a neutral resonance decaying to T(2S')7r'^, Z°(10610), has been obtained in 
a Dalitz plot analysis of T(5S') — )■ T(2S')7r°7r° decay. The statistical significance of the 
Z°(10610) signal is 5.3a (4.9cr including model and systematic uncertainties). Its measured 
mass, M (Z°(10610)) = 10609tg±6 MeV/c^, is consistent with the mass of the corresponding 
charged state, the Z^(10610). The Z°(10650) signal is not significant in either T(lS')7r°7r° 
or T{2S)TT^'n'^. Our data do not contradict the existence of Z^(10650), but the available 
statistics is not enough for the observation of this state. 
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